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Abstract 

We investigate numerically the transition between static equilibrium and dynamic surface flow 
of a 2D cohesionless granular system driven by a continuous gravity loading. This transition is 
characterized by intermittent local dynamic rearrangements and can be described by an order pa- 
rameter defined as the density of critical contacts, e.g. contacts where the friction is fully mobilized. 
Analysis of the spatial correlations of critical contacts shows the occurence of "fluidized" clusters 
which exhibit a power-law divergence in size at the approach of the stability limit. The results 
are compatible with recent models that describe the granular system during the static/dynamic 
transition multi-phase system. 

PACS numbers: PACS numbers: 81.35, 46.10, 05.60 
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Avalanches and debris flows are of special interest both for industrial and natural pro- 
cesses. Although many advances have been made in the physical understanding of granular 
flow both from a microscopic and a continuum point of view (see for example [|I], ^ ^, ^, ^ ) , 
the mechanisms leading to the transition of a granular material from a static equilibrium to 
a dynamic state are still unclear, and difficult to analyze experimentally. Recent hydrody- 
namic model based on phase transition theory describes the granular material as a two-phase 
material with a "solid" (static) and a "liquid" (flowing) phase 0. This transition, also advo- 
cated in others conditions [§, |], is important for the understanding of avalanche instability 
in the context of risk assessment. 

The purpose of this paper is to investigate numerically this static to dynamic transition 
in a minimal model configuration, e.g. a 2D cohesionless granular medium driven by contin- 
uously loading body force, here a continuous tilt under gravity. We find that the evolution is 
chacterized by intermittent local instabilities and can be described by the fraction of contacts 
where friction force is completely mobilized. A more detailed description involves "fiuidized" 
areas which percolate when the system reaches the stability limit. These results would be 
in favor of a description in terms of a multiphase system undergoing a phase transition. 

The contact friction is described by the classical Coulomb's law. Relative slip between two 
particles in contact can occur only when the friction force is fully activitated, i.e ft = ^fJ'fn, 
where ft and /„ are respectively the tangential and normal forces at the contact and fi is 
the contact friction coefficient. This defines the Coulomb threshold. Otherwise, no slip can 
occur and ft G [— /i/„, The numerical simulations were performed using the Contact 
Dynamics method [jlO|, [Tl|, based on a fully implicit resolution of the contact forces. This 
allows for an accurate determination of the ratio ft/fn independently of any numerical reg- 
ularization parameter. In all the simulations, the contact friction coefficient is = 0.5 and 
collisions are perfectly inelastic. 

The initial configuration of the pile is generated by random deposition of Np = 4000 disks 
in a rectangular box; the disks have a uniform distribution of diameters within the range 
[Dmin, Dmax]- Wc uscd the ratio Dmax/Dmin = 1-5, but we checked that for a polydisper- 
sity as large as Dmax/Dmin = 10 the results remain essentially unchanged. The granular 
samples prepared by random deposition have a rectangular shape with a nearly fiat surface, 
a thickness of SOD and a width of 120D, where D is the mean disk diameter. The initial 
coordination number, i.e. the mean number of contacts of a particle, is z ^ 3.6. It shows 
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only weak fluctuations (up to 5%) in the course of tilting. The results presented below were 
obtained from 15 independent runs in which the random processing of grain sizes is the only 
source of noise. 

The granular bed is slowly tilted with a constant rotation rate (0.001° per time step). 
The slope increases from 6 = 0° to the maximum angle of repose 6^ — 20° at which a surface 
avalanche occurs. Within the picture of ideal Coulomb's material, 6c is related to the in- 
ternal coefficient of friction of the pile fieff through the relation fieff = tan{9c) — 0.35 [jl2|. 
However, the evolution of the pile towards Of. is not monotonous. Indeed, we observe the 
occurence of local instabilities owing to the mobilization of friction between particles. Ini- 
tially ft at the contacts is only partially activated, namely \ft\ < fJ^fn- But upon tilting, 
a number of contacts reach the Coulomb threshold, i.e. \ft\ = fJ'fn- These contacts can 
not sustain further shear force increment and can lead eventually to a slip instability or the 
disappearance of the contact. We call them critical contacts. We describe the evolution of 
the pile in terms of the fraction u of critical contacts over a volume V and as a function of 
9: 

(1) 

where Nc and are respectively the number of critical contacts and the total number of 
contacts in V at slope angle 9. Since the density of contacts remains almost constant due 
to close-packing (A^ oc l^), z/ also represents the density of critical contacts. It characterizes 
the plastic state of the pile. 

The evolution of i'{9, VpUe) is reproducible from run to run even though there are rapid 
fluctuations of u within the system (see Fig |l|(a)). When averaging the u evolution over 
15 independant realizations (inset graph in Fig|l|(a)), we observe a regular increase of i'{9) 
from an intial value z/ = to a maximum value u ~ 0.08 at the maximum angle of repose 
9c. This indicates that a partial plastification occurs well before the stability limit of the 
pile. Such a transition is however characterized by rapid fluctuations of u. They are the 
signature of intermittent instabilities during which loss of critical contacts occurs through 
local dynamical rearrangements, as can be seen in the fluctuations of the mean kinetic en- 
ergy in Fig |l](b). This suggests that even though the density of critical contacts increases in 
the mean during the mobilization stage of the system, the population of critical contacts is 
renewed by these intermittent instabilities and a single critical contact is only of short-life 



3 



time due to its metastable state. During the intermittent evolution, the system can explore 
extreme i/-states as seen in Fig |l](a). When averaging over small 9 size windows the extremal 
states at each 9 and for all independant realizations, their evolution exhibits a well defined 
limit envelope (see Fig |I](c)). It shows an exponential convergence towards an asymptotic 
limit characterized by a critical state z/^ ~ 0.08, corresponding to the stability limit. The 
granular system can therefore be described in terms of an order parameter $ = z//z/c which 
varies from zero to one and characterizes the partially fluidized transition. 

Owing to the intrinsic geometrical and stress fluctuations of each individual realization, 
each pile can explore rare extremal states beyond the critical limit Vc in the course of its 
evolution. These states are metastable and followed eventually by local dynamical rear- 
rangements, the size of which can be characterized by the fraction AN/N of contacts lost in 
that event. When analysing the size of the local rearrangements as a function of the z/-state 
explored by the system prior to the event (see Fig |^), we found that the pile evolution is 
characterized by frequent weak rearrangements, less than 1%, that corresponds to //-states 
of the system below Uc, and rare and strong rearrangements corresponding to extremal states 
u > Uc- The intermittent rearrangements during the mobilization occurs both in the bulk 
and close to the free surface. 

The evolution of u may be also interpreted in terms of the mean separation distance A(^) 
between two critical contacts: 

where z{9) is the coordination number. During the mobilization, A decreases from +oo when 
u = to Xc = 2.5D when u = Vc- For A < Ac, i.e. v > Uc, the packing undergoes large-scale 
instabilities due to spatial correlations among critical contacts (see Fig. 

To investigate the distribution of critical contacts, we introduce the probability density 
function (pdf) -P(z/) of flnding a local state u within a contact neighbourghood. The size of 
this neighbourghood must not be smaller than the minimum length separating two critical 
contacts, of the order of D, nor larger than the typical size of spatial heterogeneities of u 
(to be analyzed later). In particular, we introduce the pdf -Pc(^) at critical contacts and the 
pdf P„c(^) at non-critical contacts. The distributions Pc and P„,c time-averaged over all the 
loading phase from 6* = to 6'^ are displayed in Fig ||. The range of local states u is quite 
wide for both distributions: Pc has a well-deflned peak at z/ ~ whereas Pnc is a decreasing 
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function with a high peak at z/ = 0. This indicates that critical contacts tend to appear 
preferentially in the vicinity of other critical contacts, and are strongly correlated spatially. 
It is worth to note here that the most probable local state is found to be the critical state 
that characterizes the whole pile at the stability limit. 

We study the spatial correlations by defining the mean state z/ in a neighbourghood of 
critical contacts as a function of the neighbourghood's size r and of the inclination angle 9 
of the pile. In practice, we explore r G [1-D, 2QD]. Denoting z/j the local state in the vicinity 
of a contact i, u{0, r) can be expressed as follows: 
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where Nc is the number of critical contacts in the pile. The evolution of u{6,r) is displayed 
in Fig. ^(a) for three values of 6. Independantly of the inclination 6, v decreases rapidly as 
a function of r as long as r ^ lOD. Beyond this distance, v remains almost constant and 
equal to the state of the pile z/(6', V^j/e)- Such an evolution of u{9,r) can be approximated 
by (see Fig|(b)): 

{z/oo(^)(l + A(^)r"(^)) ifr<10D, 
(4) 
Uoo{9) ifr>10D. 

where Voo{9) is the state of the pile at 9, and A{9) and a{9) are affine functions of 9. The 
distance L ~ lOD arises here as the characteristic length of the spatial correlations between 
critical contacts. 

The spatial correlations can be illustrated when looking at two snapshots of the pile for 
different values of the inclination angle 9 (see Fig In particular, the system exhibits 
clusters of critical contacts where v > Uc- These clusters can no longer sustain shear load 
increment and eventually lead to local rearrangements. They can be analysed as "fiuidized" 
zones with 9 increasing size. 

We define Tc the mean size of the clusters of critical contacts characterized by a state 
u = Uc- Equation (H) yields 



rci9) 



— - 1 



l/a{e) 



(5) 



Ai9) \u, 

where is defined as the state of the pile z/(6', VpHe) averaged over all the simulations (see 
inset graph in Fig p. The evolution of rc{9), displayed in Fig H, shows a slow increase from 



around 2D to AD as long as ^ 15 degrees. For larger 6, the size increases quickly with a 
power-law divergence oc {6c — 6)^^. That is reminescent of a percolation-like process. The 
value 9 = 15 degrees may be related to the dynamical angle of repose of the pile, but this 
would require further investigation. 

The mobilization of the granular pile is related to the occurence of "fluidized" clusters of 
critical contacts. The size of these clusters increases with 6, with a power-law divergence as 
6 6c. This result is in support of a description of the granular system during mobilization 
as a multi-phase system. The transition from a static equilibrium to a dynamical flow is 
controlled by the order parameter $ = u/uc, and the stability threshold would result from 
a multi-phase instability. Detailed analysis of the phases interaction during transition still 
requires further analysis. 

Forces in granular media have been shown to obey a broad distribution and to define a 
strong and a weak contact network. Forces transmitted along the strong contact network 
exceed the mean force in the media and are responsible for the mechanical strength of the 
medium, with a typical correlation length C, — WD [1^. By contrast, forces transmitted 
along the weak contact network mainly contributes to an average pressure. In an extended 
forthcoming paper, we show that critical contacts appear mostly within the weak contact 
network. They correspond to lower pressure as a result of a screening effect by the strong 
contact network that support the deviatoric load and controls the static stress evolution. 
Critical contacts tend to be organized as "fluidized" clusters that eventually destabilize 
the pile when their size becomes comparable with the correlation length of strong forces 
^. The avalanche instability appears as a two phase instability triggered by the interaction 
between these correlated clusters and a solid skeleton made of the strong contact network. 
The influence of material parameters such as the coefficient of friction and the degree of 
configurational disorder in the pile are the purpose of further work. The extension of the 
problem to three dimensions has to be investigated. 

We gratefully acknowledge S. Roux and D. Sornette for useful discussions and suggestions 
about this work. This research was supported by the French Ministry of Research and the 
CNRS through the ACI "Catastrophes Naturelles" and the GdR "Milieux Divises" . 
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FIG. 1: (a) Density v of critical contacts in the pile as a function of the tilt angle for a single run 
(main curve) and averaged over 15 independant runs (inset curve), (b) Mean translational kinetic 
energy {kg.m? .s^"^) of particles as a function of 9 for the same single run, (c) mean envelope of 
the maximum values of u evaluated over the 15 runs. 

FIG. 2: Proportion AN/N of contacts lost in the packing during rearrangements as a function of 
f. 

FIG. 3: Probability density function (pdf) of local states u in the neighbourghood of critical 
contacts (filled circles) and non critical contacts (opaque circles), for a neighbourghood radius 
r = 6D. 
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FIG. 4: Mean state v in the neighbourhood of a critical contact as a fonction of the size r of 
the neighbourhood (a) in hnear scale for three values of the tilt angle = 6° (solid line), = 8° 
(dashed line) and B = 12° (dotted line), and (b) in logarithmic scale for = 6° (solid line) shown 
with a line with a slope a = —1.5 (dot-dashed line). 

FIG. 5: Density of critical contacts in the pile at ~ 5° (upper picture) and at 2± 16° (bottom 
picture). Areas where v>_Vf. are in black, whereas the state i/ = are in white. 

FIG. 6: Evolution of the mean size Tc of clusters of critical contacts (see text) as a function of the 
tilt angle in linear scale (main curve) and in logarithmic scale (inset curve), where 0c = 19.96 
degrees and the line is drawn with a slope — /3 ~ 0.5. 
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